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We search for CP violation in a sample of 4.7 x 10 4 Cabibbo suppressed D° —>■ K + K~tt + tt~ 
decays. We use 470 fb^ 1 of data recorded by the BABAR detector at the PEP-II asymmetric-energy 
e + e~ storage rings running at center-of-mass energies near 10.6 GeV. CP violation is searched for in 
the difference between the T-odd asymmetries, obtained using triple product correlations, measured 
for D° and D° decays. The measured CP violation parameter is At = (l-0±5.1 s t a t±4.4 svs t) x 10 -3 . 

PACS numbers: 13.25.Ft, 11.30.Er 



In the Standard Model (SM) of particle physics, CP 
violation arises from a complex phase in the Cabibbo- 
Kobayashi-Maskawa quark mixing matrix Physics 
beyond the SM, often referred to as New Physics (NP), 
can manifest itself through the production of new parti- 
cles, probably at high mass, or through rare processes not 
consistent with SM origins. SM predictions for CP asym- 
metries in charm meson decays are generally of O(10 -3 ), 
at least one order of magnitude lower than current exper- 
imental limits 0]. Thus, the observation of CP violation 
with current sensitivities would be a NP signal. Among 
all hadronic D decays, singly Cabibbo suppressed decays 



are uniquely sensitive to CP violation in c — > uqq transi- 
tions, effect not expected in Cabibbo favored or doubly 
Cabibbo suppressed decays Q. 



In this paper we report a search for CP violation in the 
decay D — » K + K~it + tt~ using a kinematic triple prod- 
uct correlation of the form GY = Pi • (P2 x P3), where 
each pi is a momentum vector of one of the particles in 
the decay. The product is odd under time-reversal (T) 
and, assuming the CPT theorem, T-violation is a signal 
for CP-violation. Strong interaction dynamics can pro- 
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duce a non-zero value of the At asymmetry, 

At = n ° T >0) - T{Ct < 0) (1) 

T ~ r{c T > o) + r(c T < o) ' 1 ' 

where V is the decay rate for the process, even if the weak 
phases are zero. Defining as At the T-odd asymmetry 
measured in the CP-conjugate decay process, 

A T ^ n-CT>0)-r(-Cr<0) 

r(-c T > o) + r(-c T < o)' 

we can construct: 

At = ^(A t -At~), (3) 

which is a true T- violating signal Q . At least four parti- 
cles are required in the final state so that the three used 
to define the triple product are independent 0] of each 
other. Singly Cabibbo suppressed decays having rela- 
tively high branching fractions and four different parti- 
cles in the final state, therefore suitable for this type of 
analysis, are D° — > K + K~ir + n~ (explored in this paper) 
and D+ — > K+K®7r+7r~ . A full angular analysis of these 
D decays is suggested as a method for searching for CP 
violation 

Following the suggestion by I.I. Bigi 7] to study CP 
violation using this technique, the FOCUS collabora- 
tion made the first measurements using approximately 
800 events and reported At(D° -t K+K~ 
0.010 ± 0.057 ± 0.037 0]. We perform a similar study 
using approximately 4.7 x 10 4 events. 

This analysis is based on a 470 fb _1 data sample 
recorded at the T(AS) resonance and 40 MeV below 
the resonance by the BaBAR detector at the PEP-II 
asymmetric-energy e+e~ storage rings. The BABAR de- 
tector is described in detail elsewhere We mention 
here only the parts of the detector which are used in 
the present analysis. Charged particles are detected and 
their momenta measured with a combination of a cylin- 
drical drift chamber (DCH) and a silicon vertex tracker 
(SVT), both operating within the 1.5 T magnetic field 
of a superconducting solenoid. The information from a 
ring-imaging Cherenkov detector combined with energy- 
loss measurements in the SVT and DCH provide identi- 
fication of charged kaon and pion candidates. 

The reaction [lCj 

e+e- ->■ X D*+; D*+ -> n+D°; D° -> K+R-tt+tt-, 

(4) 

where X indicates any system composed by charged and 
neutral particles, has been reconstructed from the sam- 
ple of events having at least five charged tracks. We first 
reconstruct the D° candidate. All R+R~tt+7t~ com- 
binations assembled from well-measured and positively 
identified kaons and pions are constrained to a common 
vertex requiring a % 2 fit probability greater than 0.1%. 



To reconstruct the D*+ candidate, we perform a vertex 
fit of the D° candidates with all combinations of charged 
tracks having a laboratory momentum below 0.65 GeV/c 
(tt+) with the constraint that the new vertex is located 
in the interaction region. We require the fit probability 
to be greater than 0.1%. 

We require the D° to have a center-of-mass momen- 
tum greater than 2.5 GeV/c. This requirement removes 
any D° coming from B decays. We observe a contamina- 
tion of the signal sample from D a — > R+ K~ K®, where 
K® —> tt+tt" . The tt+tt" effective mass shows, in fact, 
a distinct mass peak, which can be represented by a 
Gaussian distribution with a = 4.20 ± 0.26 MeV/c 2 , and 
which accounts for 5.2% of the selected data sample. We 
veto Kg candidates within a window of 2.5 a. This cut, 
while reducing to negligible level the background from 
D° — » K + K~ K®, removes 5.8% of the signal events. 

We look for backgrounds from charm decay modes 
with mis-identified pions by assigning alternatively the 
pion mass to both kaons. Then we study the two- 
body, three-body, four-body and five-body mass distri- 
butions (including the We observe a signal of 
D+ —> R+R^tt+tt^tt+ in the five-particle mass distri- 
bution, which is taken into account in the following fit. 
No other signal is observed in the resulting mass spectra. 

We define the mass difference Am as: 

Am = m(K+K- 7r + n-7r+) - m(K+ K~ n+ 'tT ") . (5) 

Figure [lja) shows the scatter plot m(K+ K~ tx+ n^) 
vs. Am for all the events. Figure [U(b) shows the 
m(K+ K^tt+tt-) projection, Fig.QJc) shows the Am pro- 
jection. 

We perform a fit to the m(K+ K~ tx+ 7i~) and Am dis- 
tributions, using a polynomial background and a single 
Gaussian. The fit gives a D o = 3.94 ±0.05 MeV/c 2 for the 
D° mass and a D ,+ = 244±20keV/c 2 for the Am. We de- 
fine the signal region within ±2ajja and ±3.5<r D »+. The 
total yield of tagged D° mesons in the signal region is 
approximately 4.7 x 10 4 events. 

The D° yields to be used in the calculation of the 
T asymmetry are determined using a binned, extended 
maximum-likelihood fit to the 2-D (m(K+ K^tt+tt^), 
Am) distribution obtained with the two observables 
m{K+ K~7r+7T~) and Am in the mass regions defined 
in the ranges 1.825 < m(K+ R-tt+tt-) < 1.915 GeV/c 2 
and 0.1395 < Am < 0.1545 GeV/c 2 respectively. Events 
having more than one slow pion candidate in this mass 
region are removed (1.8 % of the final sample). The final 
2-D distribution contains approximately 1.5 x 10 5 events 
and is divided into a 100 x 100 grid. 

The 2-D (m(K+K-7r+n-), Am) distribution is de- 
scribed by five components: 

1. True D signal originating from a D*+ decay. This 
component has characteristic peaks in both observ- 
ables m(K+ K~7T+7r~) and Am. 
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FIG. 1: (a) m(K + K~tv + -k~) vs. Am for the total data sample, (b) m{K + K~ -k + -k~) and (c) Am projections with curves from 
the fit results. Shaded areas indicate the different contributions. The fit residuals, represented by the pulls, are also shown 
under each distribution. 



2. Random 7r+ events where a true D° is associated 
to an incorrect 7r+ , called D° peaking. This contri- 
bution has the same shape in m(K + K~ir + ir~) as 
signal events, but does not peak in Am. 

3. Misreconstructed D° decays where one or more of 
the D° decay products are either not reconstructed 
or reconstructed with the wrong particle hypothe- 
sis, called Am peaking. Some of these events show 
a peak in Am, but not in m(K + K~tt + tt~). 

4. Combinatorial background where the K + , K~, 
7r + , 7T _ candidates are not fragments of the same 
D° decay, called combinatoric. This contribu- 
tion does not exhibit any peaking structure in 
m(K + K~ tt + n~ ) or Am. 

5. Df —> K + K~TT + ir~ir + contamination, called Df. 
This background has been studied on Monte Carlo 
(MC) simulations and shows a characteristic lin- 
ear narrow shape in the 2-D (m(K + K~it + it~), 
Am) distribution, too small to be directly visible 
in Fig. QJa). 

The functional forms of the probability density func- 
tions (PDFs) for the signal and background components 
are based on studies of MC samples. These events 
are generated using the Geant4 program llj and are 
processed through the same reconstruction and analysis 
chain as the real events. However, all parameters related 
to these functions are determined from two-dimensional 
likelihood fits to data over the full m(K + K~tt + it~) vs. 
Am region. We make use of combinations of Gaussian 
and Johnson SU [T3| lineshapes for peaking distributions, 
and we use polynomials and threshold functions for the 
non-peaking backgrounds. 



The event yields and fractions of the different compo- 
nents arising from the fit are given in Table U and shown 
in Fig. [T] The fit residuals shown under each distribution 
are represented by Pull = (N da ta - N fit) I 'VNdata- 

TABLE I: Fitted number of events for each category. 



Category 


Events 


Fraction (%) 


1. Signal 


46691 ± 241 


30.8 ±0.3 


2. D° peaking 


5178 ± 331 


3.4 ±0.2 


3. Am peaking 


57099 ± 797 


37.7 ±0.6 


4. Combinatoric 


40512 ± 818 


26.7 ±0.6 


5. D+ 


2023 ± 156 


1.3 ±0.1 


Total 


151503 ± 1223 





Using momenta of the decay particles calculated in the 
D rest frame, we define the triple product correlations 
Ct and Ct as 



Ct =Pk- ■ ifw- x_Ptt+)- 



(G) 



According to the D* + tag and the Ct variable, we 
divide the total data sample into four subsamples, defined 
in Table [TTJ These four data samples are fit with fixed 
PDFs from the total sample. The signal event yields are 
given in Table ILT1 Fig. [5] shows the K + K~ ir + ir~ mass 
distributions for the four different Ct subsamples with fit 
projections in the Am signal region previously defined. 

We validate the method using e + e~ — > cc MC simu- 
lations, where D° decays through the intermediate res- 
onances with the branching fractions reported in the 
PDG [13}. We obtain a T asymmetry A T = (2.3 ±3.3) x 
1CP 3 , consistent with the generated value of 1.0 x 1CP 3 . 
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TABLE II: Definition of the four subsamples and the event 
yields from the fit. 



asymmetries: 



Subsample 


Events 


(a) D u , C T > 

(b) D°, C T < 

(c) D°, C T > 

(d) D°, C T < 


10974 ± 117 
12587 ± 125 
10749 ± 116 
12380 ± 124 



o 
> 



1000 
800 
600 
400 
200 



0- 



sou 



600 



■ data 
■□Signal 
j Background 



D°(C t >0) 





D° (C t <0) 



c 

% 

Hi 



1.84 1.86 1.88 1.9 

m(K + KV7i-) (GeV/c 2 ) 


1.84 1.86 1.88 1.9 

m(K*K'jtV) (GeV/c 2 ) 




■ D° (C T <0) ■ 














1.84 1.86 1.88 1.9 

m(K + KVjr) (GeV/c 2 ) 

: Fit projections onto the m 


1.84 1.86 1.88 1.9 

m(K + K Vtt) (GeV/c 2 ) 

[K+K--K+H-) for the four 



FIG. 

different Ct subsamples with cut on Am. The shaded areas 
indicate the total backgrounds. The fit residuals, represented 
by the pulls are also shown under each distribution. 



To test the effect of possible asymmetries generated 
by the detector, we use signal MC in which the D° 
decays uniformly over phase space. In this case possi- 
ble asymmetries are generated only by the detector ef- 
ficiency. These reconstructed events give an asymmetry 
At = (1.1 ± 1-1) x 10~ 3 , again consistent with zero. 

To avoid potential bias, all event selection criteria are 
determined before separating the data into the four sub- 
samples of Table HU Systematic uncertainties are ob- 
tained directly from the data. In these studies the true 
At and At central values are masked by adding unknown 
random offsets. 

After removing the offsets, we measure the following 



A T = (-68.5±7.3 sta1 : r,.s 
A T = (-70.5 ± 7.3 S 



Wst) x 10 ~ 3 ' 
!stat±3.9 syst )xl(r 3 . 



(7) 



We observe non-zero values of At and At indicat- 
ing that final state interaction effects arc significant in 
this D° decay. No effect is found, on the other hand, 
in the analysis of MC samples. Final state interactions 
effects are common in hadronic D decays because of the 
complex interference patterns between intermediate res- 
onances formed between hadrons in the final states [14| . 

The result for the CP violation parameter, At, is 



^ r = (1.0±5.1 stat ±4.4 8yst ) x 10- 



(8) 



The sources of systematic uncertainties considered in 
this analysis are listed in Table IIII1 The estimates of 
their values are derived as follows: 

1. The PDFs used to describe the signal are modified, 
replacing the Johnson SU function by a Crystal Ball 
function [151 ] , obtaining fits of similar quality. 

2. As in 1., for the peaking background. 

3. We increase the number of bins of the 2-D 
(m(K + K~n + iT~), Am) distribution to a (120 x 
120) grid and decrease to a grid of (80 x 80). 

4. The particle identification algorithms used to iden- 
tify kaons and pions are modified to more strin- 
gent conditions in different combinations. We no- 
tice that the difference between different selection 
efficiencies is significantly larger than the uncer- 
tainties on efficiency of the default selection. On 
the other hand, the use of the discrepancy between 
data and MC obtained using high statistics control 
samples, gives a much lower contribution. 

5. The p*(D°) cut is increased to 2.6GeV/c and 
2.7GeV/c. 

6. We study possible intrinsic asymmetries due to the 
interference between the electromagnetic e + e~ — > 
7* — > cc and weak neutral current e + e~ — > Z Q — > 
cc amplitudes. This interference produces a D° /D° 
production asymmetry that varies linearly with the 
quark production angle with respect to the e~ 
direction. Since BABAR is an asymmetric detec- 
tor, the final yields of D° and D° are not equal. 
We constrain the possible systematics by measur- 
ing At in three regions of the center-of-mass D° 
production angle 6*: forward (0.3 < cos(f5*) D o), 
central (—0.3 < cos(f5*) D o < 0.3), and backward 
(cos(6*) D o < -0.3). We observe that the At an- 
gular variation is, within the large statistical errors, 
consistent with zero as expected from the MC 
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7. Fit bias: we use MC simulations to compute the 
difference between the generated and reconstructed 
A T - 

8. Mistag: there are a few ambiguous cases with more 
than one D* in the event. We use MC simulations 
where these events are included or excluded from 
the analysis. This effect has a negligible contribu- 
tion to the systematic uncertainty. 

9. Detector asymmetry: we use the value obtained 
from the MC simulation where D° decays uniformly 
over the phase space. 

In the evaluation of the systematic uncertainties, we 
keep, for a given category, the largest deviation from 
the reference value and assume symmetric uncertainties. 
Thus, most systematic uncertainties have a statistical 
component, and are conservatively estimated. In con- 



TABLE III: Systematic uncertainty evaluation on At, At, 
and At in units of 10 -3 . 



Effect 


At 


A T 


A T 


1. Alternative signal PDF 


0.2 


0.3 


0.2 


2. Alternative misreconstructed D° 


PDF 0.5 


0.1 


0.9 


3. Bin size 


0.2 


0.4 


0.3 


4. Particle identification 


3.5 


4.2 


2.9 


5. p*(D°) cut 


1.7 


1.6 


2.4 


6. cos 6* dependence 


0.9 


0.0 


0.2 


7. Fit bias 


1.4 


3.0 


0.3 


8. Mistag 


0.0 


0.0 


0.0 


9. Detector asymmetry 


1.1 


2.1 


0.0 


Total 


4.4 


5.8 


3.9 



elusion, we search for CP violation using T-odd correla- 
tions in a high statistics sample of Cabibbo suppressed 
D° — > K + K~ir + ir~ decays. We obtain a T-violating 
asymmetry consistent with zero with a sensitivity of ~ 
0.5 %. 

The study of triple product correlations in B decays 
shows evidence for final state interaction but also give 
asymmetries consistent with zero, in agreement with SM 
expectations [3]. These results constrain the possible 
effects of New Physics in this observable Q . The results 
from this analysis fix a reference point, since the study of 
T-odd correlations play an important role in the Physics 
program of present and future charm and .B-factories. 
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